Experimental and theoretical study of processes of formation and growth of pearlite 

colonies in eutectoid steels 
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We describe our optical and electron-microscopy observations of pearlite structures in eutectoid 
steels which seem to imply that the mechanisms of formation of pearlite colonies in these steels 
differ from those observed earlier for non-eutectoid steels. A simple theoretical model to study 
kinetics of pearlite transformations is suggested. Simulations of growth of pearlite colonies based 
on this model reveal that for the volume carbon diffusion mechanism usually-supposed such growth 
is always unstable, and the steady-state growth can be realized only via the interfacial carbon 
diffusion mechanism. A model of formation of pearlite colonies based on the assumption of a 
strong enhancement of carbon diffusion near grain boundaries is also suggested. The model can 
be applicable to the plastically deformed steels, and the results of simulations based on this model 
qualitatively agree with some microstructural features of formation of pearlite colonies observed in 
such steels. 

PACS numbers: 05.70.Fh; 05.10.Gg 



I. INTRODUCTION 

Studies of growth and formation of pearlite colonies in 
steels attract great attention for many decades, see, e. g. 
1]- 10] and references therein. However, both the exper- 
imental information and the theoretical understanding 
of these phenomena seem to remain to be rather lim- 
ited. Most of experimental data used in the literature 
(see e. g. [H, H|, Q-[l(|) were obtained long ago and 
by not modern methods, and many important phenom- 
ena, in particular, formation of pearlite colonies, seem 
to be insufficiently studied. In theoretical treatments, 
only steady-state growth of colonies is usually considered, 
basing mainly on the equilibrium thermodynamics ideas 
but with no attempts of microscopic treatments 
of kinetic processes. Recently, phenomenological phase- 
field approaches were used to treat the pearlite growth 
problem [H,[l(|. However, these treatments employ many 
phenomenological parameters whose physical meaning is 
not always clear (while no microscopic estimates for them 
are given) , and they discuss only volume but not interfa- 
cial diffusion mechanisms (which can be very important 
for such processes, see [ll| and below). Thus the results 
obtained can hardly significantly elucidate the actual mi- 
croscopic mechanisms of growth of pearlite colonies. 

At the same time, even main mechanisms of pearlite 
transformations seem to be not well understood as yet. It 
is unknown whether the colony growth kinetics is deter- 
mined by the volume or the interfacial diffusion of carbon 
atoms, and this problem is debated @, [Hj|. ^ * s n °t cl ear 
which kinetic or thermodynamic factors determine the 
colony period S [l|, Q . There are absent not only the- 
ories but even any definite ideas about the mechanism 
of formation of colonies near grain boundaries and other 
lattice defects [l|, 0, Hjj, etc. 

This work aims to study some of the problems men- 



tioned both experimentally and theoretically. In sec. 2 
we present some results of our experimental observations 
of processes of formation of pearlite colonies in an eutec- 
toid steel. These processes are widely discussed in the 
literature. It is generally believed that the colonies are 
formed near grain boundaries of austenite or near other 
defects of crystal lattice [l|, , but the detailed experi- 
mental information about these processes is rather scarce 
as yet. Most thoroughly they have been discussed by 
Hillcrt who studied formation of pearlite colonies in 
both hypoeutectoid steels (those with the carbon concen- 
tration c lower than its eutectoid value c e ) and hypereu- 
tectoid steels with c > c e . Hillert found that the earlier- 
supposed "repeated sidewise growth" mechanism [l[ is 
not observed in his experiments; instead, "the individ- 
ual lamellae form by branching during edgewise growth" 
near various lattice defects: grain boundaries, twinned 
planes, interphase boundaries, etc. However, the eutec- 
toid steels with c ~ c e were not studied in Ref. 0. 
Therefore, it seems interesting to study the colony for- 
mation processes in such steels by the modern methods. 
In Sec. 2 we describe such studies and compare the re- 
sults of our observations to those of previous work [l|, Q ■ 

The rest part of this work describes some attempts to 
develop the microscopic approach to the pearlite trans- 
formation theory. As a possible first step for that, a sim- 
ple model of alloys iron-carbon was recently suggested 
hj. This is a binary interstitial alloy model with the 
phase diagrams illustrated by Fig. 5(a) below being 
symmetric with respect to the interchange of ferrite and 
cemcntite. The pearlite growth kinetics for this model 
was studied basing on the appropriate Ginzburg-Landau 
functional. The volume diffusion of carbon was supposed 
to be the main kinetic mechanism of this growth, in ac- 
cordance with the most of models used in the literature 
l]-[5]. The results of simulations of growth of colonies 
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made in [13] agreed basically with the conclusions of phe- 
nomenological treatments [3j]-[5| , though some new effects 
have also been found. 

However, the real phase diagram ferrite-cementite (if 
we model it by a binary alloy model) is by no means 
symmetric with respect to the interchanges of ferrite and 
cementite. For example, if we take the carbon concen- 
tration in cementite for unity, then the eutectoid concen- 
tration is actually close to c e =l/8, instead of c e =l/2 
in symmetric models. Therefore, to make such approach 
more realistic, we should study the colony growth also 
for non-symmetric models with c e ~ 1/8. This is done 
in Sees. 3 and 4 of the present work. We find that for 
all models of the type considered, the growth of colonies 
via the volume carbon diffusion mechanism is unstable, 
and the steady-state growth can be realized only via the 
interfacial diffusion of carbon. It agrees with the similar 
conclusion of a recent empirical analysis Q- 

In Sec. 5 we suggest a model of formation of pearlite 
colonies near grain boundaries of austenite. As men- 
tioned, the colonies are usually formed just in this region, 
but until now there seem to be no microscopic explana- 
tions for that [2,12]. In the model suggested, we relate the 
formation of colonies to a great enhancement of carbon 
diffusion near grain boundaries discussed by a number of 
authors [13, EH ■ The results of our treatment show that 
in the simplest form presented this model can be applied 
only to the strongly deformed materials, such as those 
studied in Ref. [6]. However, further generalizations can 
extend the region of applicability of the approach sug- 
gested. Main conclusions are summarized in Sec. 6. 



II. OBSERVATIONS OF MICROSTRUCTURE 
OF PEARLITE COLONIES IN NON-DEFORMED 
AND DEFORMED EUTECTOID STEELS 

A. Preparation of samples and methods of 
experiments 

All microstructural investigations were made on sam- 
ples of steel 85 with the composition shown in Table 1. 
Thus, the concentration of carbon was close to the eu- 
tectoid one for alloys iron-carbon: cf 0_c ~ 3.45%. It 
differs the steel used in our study from the hypoeutec- 
toid steels with c < c e and hypereutectoid steels with 
c > c e investigated by Hillert [2j. 

Table 1. Composition of samples. 



Impurity 


C 


Mn 


Si 


Cr 


Ni 


Cu 


S 


P 


Al 


at. % 


3.91 


0.60 


0.48 


0.08 


0.08 


0.04 


0.04 


0.02 


0.02 



Rods of steel 85 of diameter 8 mm and length about 
500 mm were subjected to various thermal and mechani- 
cal treatments described below. These different kinds of 



treatment are numbered by the same way as figures 1-4 
in which we show the microstructure of samples for each 
kind of treatment. 

(1) Heating in a muffled furnace up to 1050°C, cooling 
at air up to 650°C, putting sample in a furnace heated 
up to 650°C, annealing at this temperature for 20 sec, 
then quenching in water. 

(2) Heating in a muffled furnace up to 1000°C, cool- 
ing in the melt of sodium nitrate NaN03 at temperature 
443°C up to reaching the melt temperature, then quench- 
ing in water. 

(3) Heating in a muffled furnace up to 1005°C, cooling 
at air up to 750°C , plastic deformation at this tempera- 
ture in a torsional type plastometer with the deformation 
rate 0.2/sec up to the strain 23% (for the surface layer), 
cooling in the melt of NaNC>3 at temperature 443°C up to 
reaching the melt temperature, then quenching in water. 

(4) Heating in a muffled furnace up to 1000°C , cooling 
in the melt of NaNC-3 at temperature 443°C up to tem- 
perature 750°C , then quenching from 750°C in water. 

Samples for microstructural studies were cut in the 
shape of thin washers parallel to the cross-section of the 
rod. Then they were subjected to a many-stage grind- 
ing and polishing procedure followed by etching in a 5% 
solution of nitrogen acid in ethyl alcohol (nital) . 

For a sample treated according to regime (1), the 
photograph shown in Fig. 1 has been obtained using 
the structure analyzer Nikon Epiphot-TME and camera 
Panasonic WV-CP610/G. 

Microstructure of samples treated according to regimes 
(2)- (4) was studied by the scanning electron microscopy 
(SEM) methods. Cuts prepared for the optical mi- 
croscopy studies were used for SEM, too. However, these 
cuts were subjected to the additional ion etching. It pro- 
vides a high cleanness of the surface and enables one to 
reveal fine details of the structure |13| . The high resolu- 
tion scanning electron microscope " SUPRA 55" by LEO 
(ZEISS) with "Gemini" field emission column was used 
to acquire the images shown in figures 2-4. The in-lens 
secondary electron detector combined with a low accel- 
erating voltage (1.5 keV) was found to provide the best 
quality images. Given that the plates of pearlite intersect 
the surface of the cut at random angles, no sample tilt 
was required. 

The SEM method was chosen because its high reso- 
lution makes it to be particularly suitable for the met- 
allographic studies of alloys with a dispersed structure, 
such as the eutectoid mixtures treated in this study. The 
characteristic feature of a topographic contrast in SEM 
is the enhanced brightness of images of sharp hills and 
ledges of a surface relief, the so-called "edge-effect" re- 
lated to the increased output of electrons from such areas 
[T3 |. This enhanced brightness is seen in figures 2-4. The 
surface of pearlite areas in these figures looks as a well- 
developed relief formed by the combined action of the 
mechanical loading under polishing and the subsequent 
etching by nital 0. At the same time, martensite (ob- 
tained from the untransformed austenite) is by about an 
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order of magnitude harder than ferrite, and it does not 
contain cementite inclusions. Therefore, martensite areas 
at the cut surface are much more flat and uniform than 
those of ferrite and pearlite, and these martensite areas 
are seen in figures 4(a) and 4(b) as a grey background. 



B. Experimental results 

It is generally accepted that "in the homogeneous 
austenite, pearlite is formed practically only on grain 
boundaries" [![. Our observations agree with this con- 
clusion. In Fig. 1 we show the network of dccarburiz- 
ing formed at the surface of a sample. This network is 
formed due to the chemical interaction of carbon con- 
tained in the steel with the atmosphere oxygen, which 
results in the formation of gaseous oxides and thus in the 
removal of carbon from the steel. In Fig. 1, the original 
grain boundaries of austenite are decorated by chains of 
tiny bright grains of ferrite. Such decoration arises be- 
cause the surface decarburizing takes place primarily on 
the grain boundaries. It is so because, first, the veloc- 
ity of grain-boundary diffusion of carbon much exceeds 
that of its diffusion within the grain and, second, because 
the grain-boundary regions are much easier to access for 
oxygen as compared to inner parts of grain. Fig. 1 also 
shows that the positions of these ferrite grains seem to 
coincide with the probable initial positions of pearlite 
colonies (seen as dark areas), thus the colonies seem to 
spread just from these chains. This observation agrees 
with the above-mentioned conclusion that formation of 
pearlite colonies takes place on the original grain bound- 
aries of austenite. 

Fig. 2 shows the microstructure of a classical tiny- 
dispersed plate-like pearlite, or sorbite, for which inter- 
plate distances, or colony periods S, are about 100-200 
nm. The well-developed relief consisting of relatively thin 
cementite plates in the ferrite matrix is formed due to the 
many-stage surface treatment (polishing and etching) de- 
scribed in Sec. 2.1 which removes the surface layer of 
ferrite but make little effect on cementite. Therefore, 
dark strips ( "cavities" ) in Fig. 2 correspond to the origi- 
nal (and removed) ferrite, while bright areas ("ridges" of 
different tilt to the cut plane), correspond to cementite. 

It seems natural to suggest that the broad approxi- 
mately vertical region of ferrite in the middle of Fig. 2, 
from one side, and the thin layer of cementite along the 
left edge of this region together with a lot of approxi- 
mately horizontal colonies adjacent to this layer, from 
other side, correspond to two opposite sides of an orig- 
inal grain boundary of austenite. Then the enrichment 
of this left region by carbon (being necessary, in partic- 
ular, to form the boundary cementite layer mentioned) 
is naturally explained by "releasing" this carbon from 
ferrite regions formed to the right of this grain bound- 
ary. Therefore, it seems most probable that the above- 
mentioned approximately horizontal colonies to the left of 
this boundary have been formed just on the boundary ce- 



FIG. 1: Network of decarburizing along original grain bound- 
aries of austenite observed for a sample prepared as described 
in point (1) of Sec. 2.1. Dark areas correspond to pearlite, 
bright inclusions within these region, to ferrite grains, and 
the grey matrix, to martensite formed after a quench of an 
untransformed austenite. 



FIG. 2: Microphotograph of a nondeformed sample prepared 
as described in point (2) of Sec. 2.1. 

mentite layer (or near it). At the same time, for pearlite 
colonies to the right of this boundary, it seems natural 
to suggest that they have been formed from the ferrite 
layer. Therefore, Fig. 2 can illustrate both two earlier- 
discussed mechanisms of formation of pearlite colonies 
P, llj], that from the primary cementite and that from 
the primary ferrite, as well as the probable close inter- 
relation of these two mechanisms in the eutectoid steel 
under consideration. 

One can also observe in Fig. 2 a number of "defect" 
lines (that is, cuts of appropriate surfaces) where differ- 
ent colonies seem to meet each other or overcome some 
obstacles. Most of these lines seem to correspond just to 
"collisions" of colonies growing from different sides of a 
grain, but some of them can also correspond to lattice 
defects, for example, to twin boundaries or subgrains of 
original austenite. 

In the plastically deformed samples, both the trans- 
formation kinetics and microstructure change consider- 
ably. The transformation rate significantly increases. In 
our experiments, after the hot treatment described in 
point (3) of Sec. 2.1, the full transformation time for a 
cylindric sample of 8 mm length decreased from 10 to 4 
seconds. The degree of order and parallelism of pearlite 
colonies are notably reduced, and many colonies are frag- 
mented, which is illustrated by frames 3(a)-3(c). How- 
ever, the main microstructural features of formation of 
colonies discussed above seem to be present in these sam- 
ples, too. For example, in frame 3(a) one can distinctly 
see the cementite layer separating ferrite and plate-like 
pearlite regions, and this layer (together with its vicin- 
ity) again seems to be the most probable place of forma- 
tion of plate-like pearlite colonies. Analogous places of 
a probable formation of plate-like pearlite colonies (be- 
ing, though, less regular) are seen in frames 3(b) and 
3(c). At the same time, on "ferrite" sides of original grain 
boundaries, the formation of plate-like pearlite colonies 
in deformed samples seems to be hampered: according 
to frames 3(a)-3(c), cementite precipitates in this regions 
have usually an irregular shape being more close to the 
globular than to the plate-like one. 



FIG. 3: Microphotograph of a deformed sample prepared as 
described in point (3) of Sec. 2.1. Frames (a), (b) and (c) 
correspond to different areas of a sample. 
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FIG. 4: Microphotograph of a sample prepared as described 
in point (4) of Sec. 2.1. Frames (a) and (b) correspond to 
different areas of a sample. 

In one of experiments with non-deformed samples, 
varying cooling and temperature conditions, we suc- 
ceeded to fix an initial stage of pearlite transformation. 
Microstructures observed in these experiments are illus- 
trated by Fig. 4 where we see some transformed regions 
surrounded by the region of martensite (grey colored) 
formed under a quench of an untransformed austenite. 
More dark areas again correspond to the original ferrite 
removed by the polishing and etching treatments men- 
tioned. In frame 4(a), in the part of the transformation 
volume adjacent to its upper boundary, we seem to ob- 
serve both formation of plate-like pearlite colonies "from 
ferrite" (via the mechanism similar to that shown in Fig. 
2 to the right of the original grain boundary), and non- 
plate-like, "globular-type" precipitates of cementite sim- 
ilar to those seen in Fig. 3. At the same time, near 
right boundaries of the transformed region we see rather 
regular plate-like colonies (tilted to the cut plane) which 
could be formed via the mechanism "from the bound- 
ary cementite layer" similar to that discussed for the left 
part of Fig. 2. In frame 4(b), analogous rather regular 
colonies are observed near lower right boundaries of the 
transformed region, thus vicinities of these boundaries 
can be the regions of formation of pearlite colonies, too. 

C. Discussion of experimental results 

The above-described observations enable us to make 
following conclusions and conjectures about mechanisms 
of formation of pearlite colonies. 

1. These observations agree with the point of view 
generally accepted [l| that in the homogeneous austenite, 
the plate-like pearlite colonies form mainly near its grain 
boundaries. 

2. At the same time, our observations seem to im- 
ply that the mechanism of this formation is complex and 
includes at least two different stages. The first stage 
seems to correspond to the decomposition of austenite 
into ferrite and cementite within some immediate vicin- 
ity of the grain boundary accompanied by a transfer of 
carbon through this boundary. On the side of this bound- 
ary enhanced by carbon, the above-mentioned cementite 
layer is usually formed, and the ferrite layer is formed on 
the opposite side of this grain boundary. 

3. Then the plate-like colonies start to form, being usu- 
ally approximately normal to the original grain bound- 
ary. The mechanisms of this formation are not clear as 
yet. However, these mechanisms seem to be certainly 
different for the "cementite" and "ferrite" sides of this 
boundary, that is, for the regions enriched and depleted 
in carbon. It seems to follow, in particular, from a quite 
different morphology of new-formed colonies in these two 



regions illustrated by Fig. 2. These differences become 
still more evident for the plastically deformed samples 
when the plate-like structure of pearlite colonies usually 
remains only on the "cementite" side of the boundary, 
while on its "ferrite" side, precipitates of cementite have 
often an approximately globular rather than the plate- 
like structure. 

4. The above-mentioned conclusions about the "many- 
step" character of formation of pearlite colonies agree 
qualitatively with the similar conclusions made by Hillert 
[2j. However, he studied the hypoeutectoid and hypereu- 
tectoid steels rather than the eutectoid ones considered 
in the present work. In addition to that, Hillert studied 
mainly processes of formation of pearlite due to the in- 
teraction of initial lamellas of ferrite or cementite with 
various lattice defects, such as twin boundaries, grain 
boundaries, cementite layers on grain boundaries, etc. 
The main mechanism of formation of colonies observed 
in his studies was "branching" of this initial lamellas. For 
the formation of colonies near grain boundaries of the ho- 
mogeneous eutectoid austenite studied in our work, such 
branchings have not been observed. Therefore, one can 
believe that these branchings are characteristic of forma- 
tion of colonies in more complex conditions considered 
by Hillert [2 while we observe more simple processes. 



III. MICROSCOPIC MODELS FOR 
DESCRIPTION OF PEARLITE 
TRANSFORMATIONS 

As mentioned, even main mechanisms of pearlite trans- 
formation seem to be not quite clear as yet. Therefore, 
in theoretical approaches to this problem it seems desir- 
able to use only simplest models which yield a reasonable 
description of phase equilibria between austenite, ferrite 
and cementite, while for the rest include as small num- 
ber of model parameters with a clear physical meaning 
as possible. In the previous work [l2j], such model has 
been suggested to study the colony growth in alloys with 
symmetrical phase diagrams of the type shown in Fig. 
5(a). Below we generalize this model to the case of more 
realistic, non-symmetrical phase diagrams. 



A. Thermodynamic model 

We consider the interstitial alloys MeX c in which the 
number of interstices (interstitial voids) for atoms X is 
equal to the number of metallic atoms Me, as in austen- 
ite and in simplified models of cementite (l4l |. Then 
the concentration x of interstitial atoms X is related 
to the average filling of interstices c by the relation: 
c = x/(l — x), but for brevity, the term "concentration" 
will be used below for the average filling of interstices, c. 
We consider the models with equilibrium phase diagrams 
temperature T - concentration c, shown in Fig. 5, in 
which phases A, D and B are analogues of ferrite, austen- 
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ite and cementite, respectively. Phase D is treated as a 
disordered solid solution of interstitial (carbon) atoms 
in the FCC lattice described by the mean-field ("regular 
solution") approximation, while the phases A and B are 
described by the order parameters 77 £ which for the 
uniform equilibrium phases have the following values: 

A: (77 = 1,C = 0), B : ( ?7 = 0,C = 1), 

D: ( V = 0,( = 0). (1) 

For the iron-carbon alloys, the parameter 77 can be con- 
sidered to be proportional to the Bain or Kurdyumov- 
Zaks deformation realizing the FCC-BCC transforma- 
tion, and £ is an analogous parameter describing the 
transition from austenite to cementite [l3| . 

The inhomogeneous alloy states under consideration 
are described using the generalized Ginzburg-Landau 
functional F [TH, [2(| supposing that the characteristic 
inhomogeneity lengths li n h in the c(r), 77(1-) and £(r) 
coordinate dependences much exceed the lattice constant 
a: 

F = f — [f(c, V, + G(Vc, V77, VC)] • (2) 

J V a 

Here v a is volume per iron atom, and /(c, 77, £) is the 
free energy of a uniform alloy per iron atom. The gradient 
term G is a bilinear form of gradients Vc, V77 and V£ 
supposed for simplicity to be isotropic: 

G = g cc Vc 2 + g vv Vri 2 + g (c V( 2 

+2g cv VcVr 1 + 2 5cC VcVC + 2< 7t , c V?7VC, (3) 

while coefficients gij in this form are supposed to be 
constant. 

For simplicity, we also assume that in our model func- 
tional ([2]) the dependence of function /(c, 77, £) on pa- 
rameters 77 and C can be described by polynomials with 
at most the fourth order. Then, for relations ([T]) for pa- 
rameters 77 and £ in equilibrium homogeneous phases to 
hold at all temperatures and concentrations, the function 
/(c, 77, C) can be taken in the following form (analogous 
to that proposed in [l6[ to model the directional solidifi- 
cation processes): 

/(c, 77, C) = <p(c) + $i(c, 77) + <f 2 (c', C) (4) 

where c' — (1 — c), and (p(c) is the free energy per one 
interstice in a disordered alloy described for simplicity by 
the mean-field approximation: 

ip(c) = T(c In c+c' In c') - V cc' /2. (5) 

Positive values Vq used below correspond to the effec- 
tive repulsion between interstitial atoms which qualita- 
tively correctly reproduces the type of carbon-carbon in- 
teractions in the real austenite [lj])- The <J>i(c, 77) and 



$2(0', C) functions are taken in the following form: 

$i(c,r?)= A lC (77 2 /2 - 77 3 /3) 

+A 1 [t iV 2 /2 - (n + l)77 3 /3 + 77 4 /4], 
<D 2 (c',C)= A 2 c'(C 2 /2-C 3 /3) 

+A 2 [r 2 C 2 /2-(r 2 + l)C 3 /3 + C 4 /4]. (6) 

Here t\ is T/2Ta where Ta is the temperature of the 
A-D phase transition at c = 0; r 2 is T/2Tb where T-q 
is the B-D phase transition temperature at c' = (i. e., 
c = 1), and Ai, A 2 , A\ and Ai are some positive energy 
parameters. 

Values of the structure parameters 77 and £ in equi- 
librium homogeneous phases are found by minimization 
of functional (j4]) with respect to 77 and £ at the given 
concentration c. The minimization with respect to 77 
yields equation df/dij = d^i/drj — with the function 
$1 from ([5]), while the minimization with respect to 77 
yields an analogous equation with the derivative of func- 
tion $ 2 : with respect to £: 

77(l-77)(r 1 -77 + Aic/A 1 ) = 0, (7) 
C(1-C)(t 2 -C + A2C / M 2 )=0. (8) 

Choice of physically acceptable solutions of these equa- 
tions at different (T, c) values will be explained for the 
equation ([7]); treatment of Eq. (JSJ is similar. In Fig. 6 
we show the 77 -dependence of the $1 function in Eq. ((5]) 
for two physically different situations: 

< (n +X x c/A x ) < 0.5, (9) 
0.5 < (n + Xic/Ax) < 1. (10) 

In both cases, the $1(77) function has minima at 77 = 
and 1 (and a maximum at a certain intermediate value 
r] = t + Xc/Ai). However, in the case ©, the right 
minimum is below the left minimum, so that phase A 
with 77 = 1 is thcrmodynamically favorable, while in the 
case (|10j) . the left minimum is below the right minimum, 
and phase D with 77 = is thermodynamically favor- 
able. The line (t% + Xic/Ai) — 0.5 separating these two 
regions in the (c, T) plane is shown as a dashed line in 
the left part of figures 5(a)-5(c). It corresponds to the 
line of phase transitions between phases A and D in the 
absence of phase separation, i. e., at a constant, "frozen" 
concentration c. Similarly, the dashed line in the right 
part of figures figures 5(a)-5(c) shows the line of phase 
transitions between phases B and D at a "frozen" con- 
centration c. 

If we take into account possible phase separation, that 
is, if we minimize the total free energy @ also with re- 
spect to the number of Me and X atoms in each phase, we 
can write total equilibrium equations for any two phases, 
1 and 2, in the well-known form: 

(df/dv)i = (df/d V ) 2 = 0; (df/OCh = (df/dCh = 0; 
(df/dc)i = (df/dc) 2 = n; fi - juci = / 2 - /jc 2 (11) 

where fi is the chemical potential of carbon with re- 
spect to iron (2~j| . 
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FIG. 5: (a) Phase diagram temperature T - concentration 
c for symmetric model 1 described by Eqs. (|12p . Phases 
A, B and D correspond to ferrite, cementite and austenite; 
solid lines are two-phase equilibrium curves. Left or right 
dashed line shows the stability limit of phase D with respect 
to transition to phase A or phase B at the fixed concentration 
c. Circle shows the (T, c) values for which simulations of 
growth of colonies have been made, (b) : Same as in (a) but 
for model 2 described by Eqs. (|13p . (c) : Same as in (a) but 
for model 3 described by Eqs. (|14l) . 

FIG. 6: Dependence of <E>i/j4i on r\ in Eq. at the tem- 
perature and concentration falling in the range determined 
by inequalities © (solid line), and inequalities (|10l) (dotted 
line). 

Simulations of phase transformations described below 
were made mainly for three thermodynamic models: for 
the symmetric model 1 from [12j | with the following values 
of parameters in Eqs. ©-([5]): 

Ai, 2 = ;4i,2 = 5T A , T A = T B , V = 0, 

and for two non-symmetric models, 2 and 3, with the 
following parameter values: 

Model 2 : Ai = 7.5T A , A x = 5T A , A 2 = UT A , 

A 2 = 9A9T A , T B = 2T A , V o = 0, g cc = 3a 2 T A , 

g m = 0.05a 2 T A , = 0.5a 2 T A , gcn,cC,vt = 0> ( 13 ) 

Model 3 : Ai = 24T A , A 1 = 8T A , 

A 2 = 23T.4, A 2 = 30T A , T B = 3.8T A 

Vo = T A , g C c,riT)X< = 2T^a 2 ; g C n,cCM = ®-> (14) 

where a is the FCC lattice constant. Equilibrium phase 
diagrams for these models are shown in Fig. 5. Eutectoid 
values of concentration and temperature, c e and T e , for 
these models are: for model 1, c e = 1/2, T e = 0.43 T A ; 
for model 2, c e = 1/3, T e — 0.54 T A ; and for model 3,: 
c e = 1/8, T e = 0.83 T A . Note that in the phase diagram 
for model 3 shown in frame 5(c), the left binodal for the 
phase equilibrium A-B coincides with the c = axis 
within accuracy of drawing: solubility limits c s (T) are 
of the order of 10 -3 . 



B. Kinetic model 

Diffusion of interstitial (carbon) atoms in the course 
of pearlite transformations is described using the quasi- 
equilibrium kinetic equation method described in Ref. 
[19| . In the case of weakly inhomogeneous states under 
consideration, this equation for local concentration c(r, t) 
takes the continuum form [20l | : 

dc/dt = -divj; j a = -a/^2D aP Vp(8F/5c). (15) 



Here a and (3 are Cartesian indices, D a p is the diffusiv- 
ity tensor, while function 5F/6c — SF/Sc(r) (having the 
meaning of the local chemical potential of carbon atoms) 
is the variational derivative of functional © with respect 
to local concentration c(r): 

SF/Sc(t) = df(c, 77, C)/dc - 2(g cc Ac + g cv Ar, + g c( A() 

where A = V 2 is the Laplace operator. In describing the 
diffusivity D we take into account not only the usual 
volume diffusivity D„ but also possible surface contribu- 
tions D s which can be important due to the enhanced 
diffusion of carbon along incoherent interphase bound- 
aries [1, G3 • Therefore, the diffusivity includes both 
volume and interfacial terms and is written as follows: 

D aP = Df + a 2 Y, e^epsiD^^Vsv + d£V 7 CV 5 C)- 

7(5 

(17) 

Here the first term describes the volume diffusion (be- 
ing, generally, different in different phases which in our 
model can be described by the dependence of on i] 
and C). Terms with and describe the surface dif- 
fusion of carbon atoms along interfaces austenite-ferrite 
and austenite-cementite, respectively. Below we discuss 
only growth of plane pearlite lamellas lying within (y, z) 
plane. For these 2D problems, Cartesian indices a, (3, 7, 5 
in the second sum in (fT7|) are x or y, while e a p = —S(j a 
is the unit antisymmetric tensor with just two non-zero 
components: e xy =1 e yx —(-l). 

Temporal evolution of structure parameters 77 and C is 
described by the time-dependent Ginzburg-Landau equa- 
tions [20l | with two phenomenological relaxation param- 
eters, 7^ and 7^: 

dij/dt — —j v SF/ S — j v [<9$i (c, ij) I di]i](r) 
-2(g cr ,Ac + g^Aij + g v (A()}; 

dC/dt = ~ lc SF/6((r) = - 7C [9$ 2 (c', ()/d( 

-2(. 9cC Ac + 5j , c A77 + 5cc AC)]. (18) 

As discussed in detail in Ref. !12|, "phonon" relax- 
ation times T p h ~ {l^ 1 ,!^ 1 ) in Eqs. (jT5J) describe re- 
laxation of crystal lattice due to the anharmonic interac- 
tions of phonons, and they are by many orders of mag- 
nitude shorter than the "diffusional" times ~ Da 2 
that describe relaxation of concentration according to 
Eq. (fT5"|) being realized by the diffusional jumps of car- 
bon atoms between interstices. Therefore, one can ex- 
pect that in the course of pearlite transformations under 
study, structure parameters 77 and £ adiabatically fast 
follow the slowly varying distribution of carbon concen- 
tration c(r,7j), minimizing the free energy © at each 
given c(r, t). In particular, if the local concentration 
c = c(r,t) obeys inequality the r](r,t) value should 
be close to unity (that is, the lattice structure should 
be close to ferrite), while if inequality (ITU1) is obeyed, 
we should have: -q(r, t) ~ (that is, the lattice structure 
should be close to austenite), and analogously for £(r, t) . 
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Deviations from these "uniform" values of rj and C (i- 
e., from their values of "zero-order" in non- uniformity) 
arise only due to the presence of last, gradient terms in 
Eqs. (|18[) which in the Ginzburg-Landau approach used 
are supposed to be small. 

In our simulations, we describe this physical picture 
as follows. The generalized diffusion equation (jT5j for 
the concentration c = c(r, t) is considered as the main 
one. In the computations, it is replaced by its finite- 
difference analog and is solved using standard iterative 
methods. However, after each step of these iterations, 
the c(r, t) values obtained are divided into two groups 
for which either inequality (|9]) or inequality (fT0|) holds. 
For points r corresponding to the first group, we put 
r)(r,t) = 1, while for r corresponding to the second 
group, we put ?y(r, t) = 0. The analogous procedure 
is made for the C( r i0 values. After that, both parame- 
ters rj and £ start to evolve with time according to the 
"phonon" equations (TTS|) (again replaced by their finite- 
difference analogs) for the time interval (t, t + At) at the 
fixed values c(r, t). In these computations, we put for 
definiteness 7^ =7^, and the interval At < 0.57" 1 was 
found to be sufficient for the full relaxation of r\ and ( to 
their "quasi-equilibrium" values r][c(r,t)} and C[ c ( r i *)]• 
Then these relaxed 77 and ( values are put into Eq. (Tl5)) 
as the initial values for the next iteration in c(r, t), and 
so on. 



IV. SIMULATIONS OF GROWTH OF 
EUTECTOID PEARLITE COLONIES 

A. Methods of simulations 

Methods of simulations of the stead y-st ate growth were 
basically the same as those in Ref. [12j. We employed 
a simulation volume V s — (L x x L y x L z ) in the cu- 
bic lattice with periodic boundary conditions along y 
and z axes. The y axis was chosen along the colony 
growth direction, while the period L z along z axis was 
taken a, so that the growth of lamellar colonies parallel 
to the yz plane was simulated. The simulation length 
L x along x axis was taken as a half of the colony period 
S: L x = S/2, and the mirror boundary conditions along 
x axis were used both at x=0 and at x—S/2. Differen- 
tial equations (fT5"j) and (fT5|) were replaced by their finite- 
difference analogs with the space-step L s and the time- 
step t s , and these equations were solved by the standard 
Runge-Kutta method. Values L s = a and t s = 10~ 3 Td 
were usually employed where Td has the meaning of a 
characteristic time of diffusional jumps on the distance 
a, while the relations of this Td to diffusivities D v or D s 
in Eq. (|I7p are indicated below. The initial distribution 
of parameters c(r), ij(r) and £(r) (illustrated by frames 
7a, 8a and 11a) was chosen close to that expected for the 
steady-state growth as described in detail in Ref. 



12] 



Our simulations were usually performed for relatively 
low temperatures: T < 0.5T e , while actually pearlite 



FIG. 7: Evolution of eutectoid colonies via the volume dif- 
fusion mechanism for the symmetric model 1 with the phase 
diagram shown in Fig. 5(a) at temperature T = 0.25Ta, the 
period S = 32a, and the following values of the reduced time 
t' = t/r d : (a) 0, (b) 50, and (c) 400. The grey level linearly 
varies with local concentration c(r) between and 1 from 
white to black. 



FIG. 8: Same as in Fig. 3 but at S = 64a. 



transformations are realized at higher T > 0.9T e . How- 
ever, the colony periods S for such "realistic" T be- 
come rather large: S > 500a, and simulations of evolu- 
tion of colonies for such S become time-consuming. At 
the same time, the main aim of this work is elucidation 
of just main mechanisms but not quantitative details of 
transformations, while qualitative manifestations of these 
mechanisms seem to be weakly sensitive to the tempera- 
ture T values. In addition to that, in studies of the most 
important problems, such as the instability of growth of 
pearlite colonies via the volume diffusion mechanism dis- 
cussed in Sec. IIVB1 we also made checking simulations 
at higher T ~ 0.8X1=, and the results did not significantly 
change. Therefore, the main results and conclusions pre- 
sented below seem to be realistic and reliable in spite of 
all simplifications in both models and simulation param- 
eters used. 



B. Growth of colonies via the volume diffusion of 
carbon mechanism 

In Figs. 7-10 we present some results of our simulations 
of growth of eutectoid colonies (i. e. growth of colonies 
into austenite with the carbon concentration c a equal to 
the eutectoid one, c e ) via the volume diffusion of carbon 
mechanism. For these simulations we suppose that the 
surface diffusivities in Eq. (TT71) are absent: D]= D^=0, 
while the volume diffusivity is isotropic: D"P=D v 5 a p 
and it has the same value D v — a 2 /rd in all three 
phases under consideration, austenite, ferrite and cemen- 
tite. Various thermodynamic models used for these sim- 
ulations differ mainly by the eutectoid concentration c e 
value. Figs. 7 and 8 correspond to the symmetrical 
model 1 with c e —l/2, and the phase diagram shown 
in Fig. 5(a). Fig. 9 corresponds to a "not strongly" 
asymmetric model 2 with c e =l/3 and the phase dia- 
gram shown in Fig. 5(b). Finally, in Fig. 10 we show 
results for the more realistic, strongly asymmetric model 
3 with c e =l/8. 

The results presented in Figs. 7 and 8 (as well as other 
results for the growth of colonies in symmetrical mod- 
els discussed in detail in Ref. [l2|]) basically agree with 
conclusions of phenomenological theories [f| Hj]. At the 
same time, Fig. 9 shows that for asymmetric models, this 
growth becomes not quite stationary: it is accompanied 
by some oscillations in the colony front structure, where 
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FIG. 9: Same as in Fig. 7 but for the weakly asymmetric 
model 2 with the phase diagram shown in Fig. 5(b) at tem- 
perature T = 0.2Ta, period S — 60a, and the following t': 
(a) 100, (b) 300, (c) 500, and (d) 1000. 



FIG. 10: Same as in Fig. 7 but for model 3 with the phase di- 
agram shown in Fig. 5(c) at temperature T = OATa, period 
S = 48a, and following t': (a) 3, (b) 10, and (c) 12. 



the cementite and ferrite lamellas leave behind each other 
in turn. However, for a "not strongly" asymmetric model 
2, this non-stationarity seems to be not very important. 

However, for strongly asymmetric models, such as our 
realistic model 3 with c e =l/8, the analogous instability 
of the growth leads to the impossibility of formation of 
a regular structure of pearlite colonies via the volume 
diffusion mechanism. It is illustrated by Fig. 10 which 
shows that in this case, already for rather short evolution 
times: t < IOt^, the ferrite lamellas start to leave behind 
the cementite lamellas and to fuse with each other. It 
leads to the isolation ( "divorcing" ) of cementite lamellas 
from austenite and thus to the locking of their further 
growth. 

The similar locking (which in literature is sometimes 
called "divorcing of pearlite" @) had been observed in 
our simulations for all models with realistic c e ~ 1/8. 
This phenomenon was found to be insensitive to varying 
thermodynamic and gradient parameters Xi, Ai, Ta,b, 
Vo and gik in Eqs. (|14[) . as well as temperature T, and 
the diffusivity D v in ferrite and cementite. 

Such instability of growth of eutectoid colonies at small 
c e ~ 1/8 seems to be related to the most general relation 
between the diffusion time td and the diffusion length Id 
for the volume diffusion mechanism: td ~ l 2 /D v . As the 
distance Id needed for carbon atoms to diffuse through 
austenite between centers of the ferrite and the cementite 
lamellas much exceeds the cementite lamella half-width 
l c : I 2 ~ 50 I 2 >(j, the diffusion time needed for carbon 
atoms released from a growing ferrite lamella to reach ce- 
mentite lamellas much exceeds times of "consumption" of 
carbon atoms from the space before front of these cemen- 
tite lamellas. An amount of this fast consumed carbon is 
much lower than that needed for a noticeable growth of 
this highly concentrated cementite lamella. Therefore, it 
grows much slower than adjacent ferrite lamellas whose 
growth is realized via diffusion of released carbon atoms 
mainly "forward" , into the non-transformed austenite, 
rather than sideways, to the narrow and low-mobile ce- 
mentite lamellas. Thus the ferrite transformation front 
leaves behind the cementite transformation front. There- 
fore, the steady-state transfer of carbon atoms between 
ferrite and cementite lamellas (supposed to occur in the 
phenomenological treatments [3j-[5j) is actually not real- 
ized. 



FIG. 11: Evolution of eutectoid colonies for model 3 via the 
interfacial diffusion of carbon mechanism described in the text 
at temperature T = OATa and different periods S. Upper 
row: S — 44a, while the values of t' = t/Td are: (a) 0, (b) 
500, (c) 1000. Middle row: 3 = 48a, while t' are: (d) 100, 
(e) 400, and (f) 700. Lower row: S = 52a, while t' are: (g) 
100, (h) 500, and (i) 1000. 



FIG. 12: Same as in Fig. [TT]but at S = 40a and the following 
t': (a) 50, (b) 100, and (c) 120. 

C. Growth of colonies via the interfacial diffusion 
of carbon mechanism 

As the austenite-pearlite interphase boundary is basi- 
cally incoherent Q , one can expect the surface diffusivity 
along this boundary to much exceed the volume diffu- 
sivity D% of carbon in austenite, similarly to the grain 
boundary diffusivity which, according to measurements 
by Bokshtein et al. [l7j], exceeds this volume diffusiv- 
ity by several orders of magnitude. We also note that 
the volume diffusivity in ferrite at temperatures under 
consideration exceeds that in austenite by 2-3 orders of 
magnitude |22} . Keeping in mind all that and employing 
also considerations of simplicity, in simulations of growth 
of colonies via the interfacial diffusion mechanism we sup- 
posed the values of effective interfacial diffusivities, pro- 
portional to quantities and in Eq. flTTl) , and the 
volume diffusivity in ferrite, to be similar, while the 
volume diffusivities in austenite and cementite, D%, and 
D%, to be negligibly small. Therefore, in these simula- 
tions we put: 

D1 = Df = 20a 2 / Td , D( = a 2 /T d , 

D a v = D% = 0, (19) 

where we also take into account that interphase bound- 
aries in our model have widths w ~(4-5)a, thus the struc- 
ture parameters gradients Vrj and V£ in Eq. (fl~7f can 
be estimated as: \Vr)\ ~ |VC| ~ (0.2-0. 25) /a. 

Some results our simulations for model (|19p are pre- 
sented in Figs. 11-14. Let us discuss these results, First, 
Fig. 11 shows that the stable steady-state growth of 
colonies via the interfacial diffusion mechanism is pos- 
sible, unlike that via the volume diffusion mechanism 
discussed in Sec. IIVBI Fig. 11 also shows that both 
the growth rate V and the stationary front shape no- 
tably vary with the colony period S which qualitatively 
agrees with conclusions of phenomenological treatments 
0-11. At the same time, our microscopic approach re- 
veals many kinetic features which are absent in these 
treatments. In particular, Figs. 11-14 show that the 



FIG. 13: Same as in Fig. [TT]but at S — 56a and the following 
t': (a) 50, (b) 250, and (c) 280. 
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FIG. 14: (A): Ratio of the steady-state colony growth rate V 
to its maximum value V ma x versus the reduced colony period 
S/ S(Vmax)- Solid line: our simulations illustrated by Fig. 11 
for which we found: V m ax~0.031a/Td. Dotted and dashed 
lines: results of phenomenological treatments 0, [H] for the 
volume and surface diffusion mechanism, respectively. (B): 
Distribution of colony periods S observed in experiments [lj . 



growth rate V depends on the period S much stronger 
than in the phenomenological treatments, and the inter- 
val of possible periods S is limited not only from below 
(by the minimal value So which, according to Zener Q , 
is determined by the balance between the volume gain 
and the surface loss of free energy under colony growth), 
but also from above, by some maximum value S ma x re- 
lated to the development at S > S max of the pearlite di- 
vorcing processes illustrated by Fig. 13; these processes 
are analogous to those shown in Fig. 10. 

The resulting dependence V(S) (shown in Fig. 14(A) 
by solid line) turns out to be much sharper than the 
analogous phenomenological dependences (shown in Fig. 
14(A) by dotted and dashed line) for both the volume 
and the surface diffusion mechanism. At the same time, 
this more sharp dependence V(S) seems to better agree 
with the distribution of periods S observed in experi- 
ments and illustrated by Fig. 14(B). 

Let us also note that our main conclusion that the 
growth of pearlite colonies is determined by the surface 
rather than volume diffusion mechanism agrees with a 
similar conclusion made by Whiting Q basing on his 
analysis of experimental data about the pearlite growth 
velocities. 



V. MODEL OF FORMATION OF PEARLITE 
COLONIES NEAR GRAIN BOUNDARIES OF 
AUSTENITE 

As mentioned, possible mechanisms of formation of 
pearlite colonies are widely discussed in the literature 
HIB Hi- However, these discussions include usually either 
detailed phenomenological description of observations of 
these complex processes 0, or just general considera- 
tions about such mechanisms [l|, H| , with no attempts of 
microscopic treatments or modeling of these processes. 
In this section we discuss a simple theoretical model of 
formation of pearlite colonies near grain boundaries of 
austenite based on the assumption of a great enhance- 
ment of carbon diffusivity in this region. We show that 
in the simplest form described below, this model can be 
applied only to strongly deformed materials but not to 
the usual materials with relaxed grain boundaries. How- 
ever, some of results described below can also be useful 
for understanding of similar processes near relaxed grain 
boundaries, in particular, those shown in Figs. 2-4 of the 
present work. 



FIG. 15: Formation of colonies for the model described 
by relations (A)-(C) in the text for the case h,ho S> a at 
Cb = c e = 0.125, the simulation length along x axis equal 
to L x =512a, and the following values of the reduced time 
t' = tD b /a 2 : (a) 0, (b) 200, (c) 300, (d) 1000. 



FIG. 16: Same as in Fig. 15 but at c b = 0.25, L x = 240a, 
and the following t': (a) 0, (b) 50, (c) 100, (d) 300. 

A. Kinetic model 

The thermodynamics of transformation will be de- 
scribed by our most realistic model 3 with parameters 
given by Eqs. f|14|) . In treatment of kinetics, it is con- 
venient to separately consider the formation of colonies 
near grain boundaries of austenite and their subsequent 
growth inside the grain. For brevity, these two stages of 
evolution will be referred to as the "formation" and the 
"growth inside the grain" stages. In simulations of the 
formation stage we made the following assumptions. 

(A) Within a layer of a width h adjacent to a plane 
grain boundary, the volume diffusivity of carbon much 
exceeds both volume and surface diffusivities within the 
grain. In our modeling, we describe it by the following 
relation generalizing Eq. (|17[) to this non-uniform case: 

0<y<h: D a p=D b 8 a p, 

h<y: D aP = D b 5 aP {l - exp[(y - h)/l]} (20) 

where parameter / characterizes the width of transition 
to the inner part of grain. 

(B) This boundary layer is enriched by carbon and has 
concentration c b > c e while the concentration within the 
grain is c e . 

(C) At the initial time t = 0, there exist a plane 
lamella of cementite (or cementite with adjacent ferrite) 
which has a length ho, width wq, and is normal to the 
grain boundary, as illustrated by frames (a) in Figs. 15- 
19. 

These assumptions qualitatively agree with the avail- 
able experimental observations and theoretical consid- 
erations. In particular, Bokshtein et al. [17] found 
the carbon diffusivity near grain boundaries in ferrite at 
T = 550° to exceed that within the grain by 3-4 orders 
of magnitude, and they observed a similar (though some- 
what lower) enhancement of diffusion near grain bound- 
aries of austenite. Authors of a recent theoretical work 
[l8j discussed acceleration of diffusion near grain bound- 
aries in strongly deformed materials; they concluded that 
the enhancement of diffusion in such materials should 
spread for the significant distances from grain bound- 
aries. The enhanced concentration of carbon and car- 
bides in vicinities of grain boundaries was noted by a 
number of authors [l|, [TtJ • For the processes of formation 
of colonies discussed above in connection with Figs. 2 
and 3, such enhancement seems to occur on the "cemen- 
tite" sides of grain boundaries, and so on. 
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FIG. 17: Same as in Fig. 15 but at c b = 0.25, L x = 240a, 
and the following t': (a) 0, (b) 50, (c) 100, and (d) 300. 



FIG. 18: Same as in Fig. 15 but at c b = 0.25, I = 4a, 
wq — 6a, h = ho = 20a, L x — 600a, and the following t': 
(a) 0, (b) 100, (c) 500, (d) 1000, (e) 1700, and (f) 2800. 

Methods of simulations of formation of colonies for the 
model (A)-(C) were basically the same as those described 
in Sec. 4. Differences concern only boundary conditions. 
For "one-dimensional" simulations shown in Figs. 15 and 
16, we put no boundary conditions along x axis, while 
for the simulations shown in Figs. 17-20, the mirror sym- 
metry with respect to both plane x — and plane y = 
was supposed. 

After a relatively fast formation of a colony near grain 
boundary discussed above, its further growth into the 
grain was supposed to occur via the interfacial diffusion 
mechanism described in Sec. 4.3. To simulate this pro- 
cess, we used the following model. 

(D) The initial carbon concentration c(r, t=0)= co(y) 
starts to gradually decrease to its volume value c e when 
the distance y from the grain boundary exceeds the 
width j/o of the layer enhanced by carbon: 

V < Vo ■ c = c b ; 

y > y : c = c e + (c b - c e ) exp [{y - yo)/l c ]. (21) 

In our simulations we used such values of parameters: 
c b = 0.25, c e = 0.125, y a = 40a, l c = 10a. 

(E) Diffusivity for this process corresponds to the in- 
terfacial diffusion mechanism described by Eqs. (|17p and 

CD2). 

(F) As the initial state for this modeling, we used 
the distribution of parameters c(r), r\ (r) and £(r) ob- 
tained in the end of simulation of formation of colonies 
for h = ho = 13a shown in Fig. 19. This distribution 
is presented in frame 19(d) for the interval of x between 
center of the second cementite lamella (nearest to the 
initial one) and center of the third ferrite lamella, with 
mirror boundary conditions with respect to both bound- 
aries of this interval along x-axis. 

B. Results of simulations of processes of formation 
of pearlite colonies 

Some results of simulations based on the above- 
described models are shown in Figs. 15-21. Let us first 
discuss the results presented in Figs. 15 and 16 which 
correspond to a limiting case of a "very thick" layer 



FIG. 19: Same as in Fig. 17 but at h = h = 15a, L x = 480a, 
and the following t': (a) 0, (b) 100, (c) 500, (d) 1000, (e) 
1700, and (f) 2300. 



FIG. 20: Same as in Fig. 17 but at h = ho = 13a, L x = 300a, 
and the following t': (a) 0, (b) 100, (c) 1000, and (d) 1600. 



FIG. 21: Same as in Fig. 17 but at h = ho = 10a, L x = 192a, 
and the following t': (a) 100, (b) 500, and (c) 1000. 

of enhanced diffusion: h, ho 3> a, when the evolution 
becomes effectively one-dimensional. Our simulations 
showed that in this case, the perfectly periodic pearlite 
structures are formed practically at any width of initial 
lamcllas wo and any concentration c b within the layer, 
as for this geometry there are no "losses" of carbon atoms 
from the phase transformation zone. For example, when 
a ferrite lamella grows into austenite, "excess" carbon 
atoms in austenite released from ferrite are accumulated 
before the front of this lamella until their local concen- 
tration exceeds the "critical" value that corresponds to 
the right dashed curve in the phase diagram of Fig. 5(c) 
(for the given temperature ,T). Then a cementite lamella 
starts to form via the "adiabatic" mechanism described 
in Sec. 3.2. Later on this lamella thickens consuming 
the carbon atoms before its front, and so on. Note that 
when concentration c b in the layer is equal to the eutec- 
toid one, c b = c e , the period of a self-organized colony 
obtained in the simulation shown in Fig. 15: S b ~ 39a, 
is very close to the minimal period of the steady-state 
growth shown in Fig. 14(A): So — 40a. 

Figs. 17-20 illustrate evolution of microstructure for 
more realistic cases of a finite width h of a layer of en- 
hanced diffusion. We found that the type of this evo- 
lution is weakly sensitive to varying parameters I, wo 
and ho mentioned in points (A) and (C), but it sharply 
depends on the width h of the layer of enhanced diffu- 
sion. In particular, at h = 20a (Fig. 17), in the pearlite 
colony formed via this mechanism survive only cemen- 
tite lamellas formed "in one", see frame 17(f). Values 
h ~ (13-15)a (Figs. 18 and 19) seem to be "optimal" for 
formation of regular colonies, but their period S sharply 
depends on the value of width h varying from S'~35a 
at h — 15a, to S'~47a at h — 13a. Finally, in too nar- 
row layers: h < 10a (Fig. 20), pearlite colonies do not 
form. Therefore, formation of regular pearlite colonies 
via the model mechanism (A)-(C) is possible only if the 
width of the layer of enhanced diffusion notably exceeds 
interatomic distances: h > 10a S> a. 

If this condition is obeyed and a colony can be formed 
near grain boundary via the mechanism (A)-(C), its fur- 
ther growth into the grain can be described by the model 
(D)-(F) of the previous section. To simulate this pro- 
cess we used the above-described formation model with 



FIG. 22: Growth of new-formed colonies into the grain for the 
model described by relations (D)-(F) at L x — 73a L y = 73a, 
and the following t' = t/r d : (a) 100, (b) 700, b (c) 1000, 
wjere r d is the same as in Sec. 4.3. 
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FIG. 23: Role of subgrains of austenite in formation of sub- 
colonies of pearlite (interrupted pearlite transformation) in a 
plastically deformed steel. Symbol "g.b." means grain bound- 
ary of austenite, "s.b." means subgrain of austenite, and dark 
regions correspond mainly to cementite (photo 34 in |(|). 

FIG. 24: Formation of lamellas from initial globulae of ce- 
mentite on a subgrain of austenite (photo 35 in [(|). 

h=13a for which the period of new-formed colonies ob- 
served in Fig. 19: S ~ 47a, is close to that correspond- 
ing to the maximum steady-state growth rate and shown 
in Fig. 14(A): S{V max ) = 48a. In Fig. 21 we present 
some results of this simulation. One sees that the colony 
steadily grows with the rate V ~ V max which seems to 
be natural for the model used with S ~ S{V max ). 

Now let us discuss a possible relation of a simple 
pearlite formation mechanism described in this section 
to reality. As mentioned, such mechanism can be ef- 
fective only if the enhanced diffusion layer is sufficiently 
wide: h > (13 — 15)a. For the usual, relaxed grain 
boundaries such widths seem to be too high. However, 
in materials subjected to an intense plastic strain, the 
enhanced diffusion regions, according to theoretical es- 
timates [ll], should significantly broaden, and the in- 
equality mentioned can be realized. In this connection, 
it can be interesting to compare our figures 18 and 19 
to some experimental data about formation of pearlite 
colonies in strongly deformed steels obtained by Tushin- 
sky et al. [(| and presented in Figs. 22 and 23. Tushinsky 
ct al. believe that these micrographs show formation of 
pearlite colonies on subgrains of austenite arisen due to 
the intense "thermoplastic hardening" of this steel. The 
morphology of some colonies seen in these figures, in par- 
ticular, those positioned to the left of the center of Fig. 
22, seems to be rather similar to that shown in Figs. 18 
and 19. 

In the usual, not deformed steels, the pearlite colonies 
seem to be formed via more complex, many-stage pro- 
cesses discussed in Sec. 2 and Ref. Q. However, one 
can believe that the peculiar features of transformations 
with strongly inhomogeneous and anisotropic diffusivity 
described above can also be manifested in these, many- 
stage paths of formation of pearlite colonies. 

VI. CONCLUSIONS 

Let us summarize the main results of this work. 
1. The microstructure of pearlite colonies in both non- 
deformed and plastically strained eutectoid steels has 



been investigated using the optical and the scanning elec- 
tron microscopy methods. The results obtained enable 
us to make some new conclusions about mechanisms of 
formation of pearlite colonies, in particular, about their 
many-step character and about differences of these mech- 
anisms for the processes of formation of colonies near 
grain boundaries in eutectoid steels and near other lat- 
tice defects in non-eutectoid steels. 

2. A simple model for theoretical studies of pearlite 
transformations is proposed. In spite of its simplicity, it 
seems to reflect the most significant features of thermo- 
dynamics and kinetics of phase transformations between 
austenite, ferrite and cementite. 

3. Simulations of growth of pearlite colonies based 
on this model showed that the volume diffusion of car- 
bon mechanism supposed in the most of existing theo- 
ries leads to the instability of the steady-state growth of 
colonies at any parameters of the model and any tem- 
peratures. A more adequate theory of growth based on 
the mechanism of interfacial diffusion of carbon is sug- 
gested. The kinetic characteristics of growth obtained in 
this theory differ notably from those obtained in the stan- 
dard phenomenological models but they seem to better 
agree with available experimental data. 

4. A model of formation of pearlite colonies based on 
the assumption of a strong enhancement of carbon dif- 
fusion near grain boundaries of austenite has been sug- 
gested. The results of simulations of processes of for- 
mation of pearlite colonies based on this model seem to 
qualitatively agree with the available data for plastically 
deformed steels. Further generalizations of this model 
can help to understand more complex many-stage pro- 
cesses observed under formation of pearlite colonies in 
non-deformed steels. 
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